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ABSTRACT 


To illustrate the way of calculus of moment, shear and axial force diagrams 
in frames with straight haunches, this article starts showing how does these 
diagrams are found in frames with constant section elements, which it’s made of two 
ways, the first one is starting from elastic ordinates and the second one using 
equilibrium. 


The way to find the diagrams it’s developed for the analysis method 3, in 
which it’s considered a compound beam with three elements, that are: left haunch; 
the central beam with constant section and the right haunch. Shape functions for the 
right haunch are obtained, since for the left haunch was shown in Aguiar and Bonilla 
(2020, 3) and for the constant section beam it was shown in Aguiar (2020). Shape 
functions constitute the Complementary Problem. 


For the solution of the Primary Problem facing a uniform distributed load, the 
fixed beam is solved for the two cases: left haunch and right haunch, in this way 
equations for shear and moment are found. In order not to distract the lector, the 
numeric development is shown in Appendix 1. The total solution is the sum of 
Primary Problem and Complementary Problem according to the sign convention 
adopted. 


Both for frames with constant section elements and for frames with 
haunches, new functions had been developed and incorporated to the computational 
system CEINCI-LAB in which are described input data for constant section and 
similarly with data for beams with haunches, in this ones it is shown with more detail 
what do this new functions do. 


As an application it is found the group of diagrams in a frame with circular 
columns and beams with haunches, applying new functions of CEINCI-LAB and 
results were compared with the ones obtained in programs such as SAP 2000 or 
ETABS. 


Keywords: Moment, shear and axial load diagrams. Concrete haunches. CEINCI- 
LAB 
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DIAGRAMAS DE FLEXION, CORTE Y CARGA AXIAL 
EN PORTICOS CON CARTELAS DE HORMIGON 


RESUMEN 


Para ilustrar la forma de calculo de los diagramas de momento, corte y carga 
axial de porticos con cartelas, se empieza el articulo presentando como se hallan 
estos diagramas en porticos con elementos de sección constante, lo cual se lo hace 
de dos formas, la primera a partir de las ordenadas de la elastica y la segunda 
mediante equilibrio. 


La forma de encontrar los diagramas se lo desarrolla para el modelo de 
analisis 3, en que se considera a una viga con cartelas compuesta por tres 
elementos, que son: la cartela izquierda; la viga central de seccion constante y la 
cartela derecha. Se obtienen las funciones de forma para el caso de la cartela 
derecha, ya que para la cartela izquierda se lo presento en Aguiar y Bonilla (2020, 
3) y para la viga de seccion constante se tiene en Aguiar (2020). Las funciones de 
forma constituyen el Problema Complementario. 


Para la solución del Problema Primario ante carga uniforme distribuida, se 
resuelve la viga empotrada-empotrada, para los dos casos: cartela a la izquierda y 
cartela a la derecha, de esta manera se hallan las ecuaciones de cálculo del cortante 
y del momento. El desarrollo numérico se indica en el Anexo 1, para no distraer la 
atención del lector. La solución total es la suma del Problema Primario y 
Complementario, teniendo en cuenta la convención de signos con que se trabaja. 


Tanto para pórticos con elementos de sección constante como para pórticos 
con cartelas, se han desarrollado nuevas funciones que han sido incorporadas al 
sistema de computación CEINCI-LAB las mismas que se describen la entrada de 
datos para el caso de sección constante y como son similares los datos para vigas 
con cartelas, en estas se indica con un poco más de detalle que hacen las nuevas 
funciones. 


Como aplicación se encuentran los diagramas en un pórtico con columnas 
circulares y vigas con cartelas, aplicando las nuevas funciones de CEINCI-LAB y 
los resultados se comparan con los que se obtienen con el programa SAP 2000 o 
con ETABS. 


Palabras Claves: Diagramas de momento, corte y carga axial. Cartelas de 
hormigón armado. CEINCI-LAB 
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1. INTRODUCCIÓN 


Se han publicado 3 artículos sobre el análisis de estructuras de hormigón 
armado con vigas con cartelas lineales; en el primero se analiza la viga con cartela, 
en la forma que se indica en la figura 1 a, (Aguiar y Bonilla, 2020, 1); en el segundo 
artículo se considera un nudo en la mitad de la viga, de tal manera que se tienen 
dos elementos uno con la cartela a la izquierda y otro con la cartela a la derecha 
(Aguiar y Bonilla, 2020, 2); en la tercera publicación se trabaja con dos nudos 
interiores en la viga, ver figura 1 c, (Aguiar y Bonilla, 2020, 3). 
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(c) 
Figura 1 Modelos de analisis; a) Sin considerar nudo en la mitad de la viga; b) 
Considerando nudo en la mitad de la viga, con lo que se tienen 2 
elementos; c) Trabajando con tres elementos en la viga. 


Para ir avanzando en el Proyecto de Investigacion: “Estructuras con 
Cartelas de Hormigón Armado o Acero” (hitps://ceincilab.wordpress.com) se 
presenta ahora el cálculo de los diagramas de momento, corte y carga axial, en el 
modelo de análisis número 3. 


Para que el lector pueda entender con facilidad, primero se encuentran 
estos diagramas para pórticos con vigas y columnas de sección constante. Luego 
se presentarán las funciones de forma para el elemento cartela a la derecha, ya que 
para el elemento cartela a la izquierda fue publicado en Aguiar y Bonilla (2020, 3) y 
se indicara el procedimiento con el cual se hallan los diagramas, que dicho sea de 
paso son útiles para el diseño de las vigas. 


2. PÓRTICOS CON ELEMENTOS DE SECCIÓN CONSTANTE 


2.1 Funciones de Forma para elementos de sección constante 


En la figura 2, con líneas entrecortadas se presenta un elemento BC, antes 
de ser deformado por la aplicación de cualquier tipo de cargas; con línea continua 
se indica al elemento deformado, el nudo B pasó a B’; el nudo C pasa a С”. Para 
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definir las coordenadas de B' se tienen las coordenadas del elemento denominadas 
U,,V;,9,, que son la componente de desplazamiento axial y transversal, y la rotación 
del nudo inicial. Para el nudo C’, las coordenadas son u,,v.,,6., con el mismo 
significado, pero en el nudo final. En resumen, las coordenadas locales del elemento 
se designan con las letras u, v, 9, seguidas de un subíndice que es 1, para el nudo 
inicial y 2 para el nudo final. La convención de signos positiva es la indicada en la 
figura 2. Aguiar (2020). 





We) 


- лу 


Š „Ж 
„ее САР 
y al 
aa 
Ч\\ „7 ow 
DN 


Figura 2 Coordenadas locales de un elemento y ordenadas de la elástica. 


Asociado al desplazamiento axial u, se tiene la función de forma 0, (x); al 
desplazamiento transversal v, se tiene la función de forma @„(х); а la rotación 0, , 
etc. Las funciones de forma para sección constante se encuentran deducidas en el 
Capítulo 3 de Aguiar (2020). 


Resta por indicar que las funciones de forma también se llaman funciones 
de interpolación, ya que conocidos u,, v,, 6, , en el nudo inicial y u,, v,, 6, en el nudo 
final, se puede hallar las ecuaciones que definen la curva indicada en la figura 2, 
llamadas u(x), v(x), con funciones de interpolación. Aguiar (2020). 


Para elementos de sección constante, estas funciones, son: 


O(c) =1-5 (1) 
ке: ae 
Øx) = 1-377 t273 (2) 
X. 2 
$0 =x(1-7) (3) 
Zo = (4) 
bs(x) = (3-25) (8) 


x 


de (x) = -—(1-=) (6) 
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Donde L es la longitud del elemento; X es la ordenada que define la posicion 
de un punto P, que se mide axialmente a partir del nudo inicial. 


2.2 Ordenadas de la Elastica 


Al aplicar el principio de superposicion en cada una de las deformadas que 
conducen al calculo de las funciones de forma se encuentran las ordenadas de la 
elastica para el problema complementario, que se indican a continuacion. 


(7) 


u(x) = u * $i G) + u; * $, (x) 
v(x) = uQ; (x) + 0,$s (x) + ®›ф»5(х) + 0,фь(х) (8) 


Para encontrar las ordenadas de la elástica, (7), (8), se debe tener presente 
el problema primario (en el cual los elementos se encuentran empotrados) y el 
problema complementario (en el que se tienen solo cargas en los nudos). La suma 
de los dos, reportan la solución total. Aguiar (2020) 


u(x) = u, * $ó, (x) + u; * ó, (x) + u, (9) 


V(X) = 0, P2(x) + 0,ó$s (x) + u; $s G) + O26 (x) + Vy (10) 


Donde u,, v, son la solución del problema primario, de la componente axial 
y transversal al eje del elemento. Se van a encontrar los diagramas de flexión, corte 
y carga axial, en pórticos planos con cargas en los nudos y con carga uniforme 
distribuida vertical en los elementos, que es el caso más frecuente que se tiene, en 
la práctica profesional. 


Para el caso de carga uniforme perpendicular al eje del elemento, u, = 0 y 
el desplazamiento vertical, en una viga de sección constante se halla con la 
siguiente expresión, que se obtiene de la resistencia de materiales, con la 
convención de signos de que el desplazamiento vertical es positivo si va hacia abajo. 


Fo 11 
= Х# +Х?12 — 2Х31, (11) 
Up = Sap + ) 





Ahora, en la figura 2 se indica la convención de signos con que se ha venido 
trabajando, se ve que (10) es positivo si va hacia arriba, para continuar con esta 
convención de signos se cambia de signo a la solución del problema primario (11), 
con lo que se tiene: 


b 
24ЕІ 





0(х) = 01ф,(х) + 01фз(х) +о,ф5(х) + 0,ф,(х) – (Х* PALAS 2X*L) 


(12) 


Se deja constancia que la ecuación 9, es valida solo para elementos de 
sección constante con carga uniforme distribuida; si se tiene cargas puntuales 
actuando sobre la viga no es válida o si se tiene una viga con cartelas tampoco es 
valida por lo que se deben hallar las funciones de forma para ese caso. 
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2.3 Cálculo de los diagramas de carga axial, momento y corte 


Para elementos de sección constante el cálculo de los diagramas de fuerza 
axial N(x); momentos M(x); Corte V(x) es directo y se lo hace con las ecuaciones 
de la resistencia de materiales: 


N(x) =EAx*u'(x) (13) 
M(x) = E | ZO 
d M(x) 
e 


Donde E es el módulo de elasticidad del material; 4 es el área de la sección 
transversal: J es el momento de inercia. 


Por lo tanto, para hallar el diagrama de carga axial, se debe derivar primero 
u(x), para encontrar el momento se debe derivar dos veces v(x) y para el cálculo 
del cortante se deriva tres veces v(x). 


Conviene derivar en las funciones de forma por separado antes que derivar 
en la ecuación 12. Con esta indicación, para hallar el diagrama de momentos se 
tiene: 


6 12x 4 6x 


Ø; (x) = a Y 05 (х) = DES 
6 12x 2 6x 
ee. ae ee ae 


Po 
M(x) = El |v, n(x) + 01 Ø3 (x) + vz Øe (x) + 03 Ø; (x) — an x? + 21? — 12x L) (14) 


Luego la ecuación 11, es la que se programa para hallar el diagrama de 


momentos. Para el diagrama de corte se debe derivar la ecuación 14, con lo que se 
obtiene. 


12 6 12 6 
JO 000) = 00) =-7ғ Osi 


Po 
24 El 





V(x) = EI |v, O4'(x) + 0, OY’ (x) + u, BE (x) + Oy OY’ (x) — (24x —12L)} (15) 


Es importante destacar que en las funciones de forma no se ha considerado 
el efecto de corte ya que la obtencion parte de resolver la siguiente ecuacion 
diferencial, para el caso de seccion constante. 


d*v 


Exl 
е 


= 0 
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e EJEMPLO 1 


La estructura de la figura 3 a, esta compuesta por columnas de 35/35 y vigas 
de 30/30; es de hormigon armado con un modulo de elasticidad E, = 


1800000 elias sobre ella gravita una carga uniforme distribuida de P, = А En la 


figura 3 b, se indican los grados de libertad con los que se analiza la estructura. Se 
pide calcular detalladamente los momentos y cortantes en cada elemento, a partir 
del vector de coordenadas generalizas q, cuyos valores se indican en la figura 3 c. 
Primero se encontro la matriz de rigidez de la estructura K, y el vector de cargas 
generalizadas Q. Luego de la solucion del sistema de ecuaciones lineales Q = K q, 
se hallo q 


2 T/m 











Ë 
со 
-0.0000494 -0.0001309 -0.0001997 
-0.0001293/" -0.0001527 / / | -0.0000646 
> — E 
Sw bw bw 
-0.0005685 -0.0003014 (0.001191 
Е 
со 
4 т 5m 


Figura 3 a) Estructura de Ejemplo 1; b) Modelo numérico de análisis; c) Grados de 
libertad. 
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e SOLUCIÓN 


Para encontrar el diagrama de momentos lo primero que se obtiene son los 
desplazamientos y giros de cada elemento en coordenadas locales, a partir del 
vector de coordenadas generalizadas q. En la tabla 1, se presentan estos valores, 
el elemento 1 es la columna izquierda; el 2 la columna central y el 3 la columna 
derecha; el 4 la viga izquierda y el 5 la viga derecha. 


Tabla 1 Coordenadas locales en los elementos que se encuentran del vector q 


DESPLAZAMIENTOS 
# ELEM TIPO 


w |» | © | ш | > 


e 
| 1 [columna] ooo | ooo | ооо 00000494 00001293 | 00005685 


Para cada elemento se tiene las coordenadas locales, que intervienen en 
las ecuaciones 11 y 12, con lo que se hallan los valores indicados en las tablas 2 y 
3, se ha calculado cada décimo de la luz. Al graficar estos valores se halla los 
diagramas de momento y de corte. 





Tabla 2 Valores de momentos en cada décimo de la longitud de cada elemento 


О ЫЕ E 





Tabla 3 Valores de cortantes en cada décimo de la longitud de cada elemento 


s| o i oooO OSOS 





Es muy importante haber desarrollado paso a paso el cálculo de los 
diagramas de flexión y de corte, para que el usuario sepa cómo se halla el diagrama 
de momentos, corte y carga axial en el computador. Solo resta graficar, pero para 
ello se tiene los programas que se indican en el siguiente apartado, destacando que 
estos diagramas se encuentran a partir de las ordenadas de la elástica. 
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2.4  — Nuevas funciones para CEINCI-LAB. Primera parte 


e  dibujo_axial(X, Y, NI, NJ, q,VC,L,seno,cose,nudt, Fm, E, ELEM, FF, fact, nelem) 


La función dibujo_axial dibuja el diagrama de carga axial del pórtico; los 
datos de entrada son: X,Y, Vectores con las coordenadas de los nudos; MNI,NJ, 
vectores con los nudos inicial y final de los elementos, respectivamente; q vector de 
coordenadas generalizadas que contiene los desplazamientos y giros de los nudos; 
VC matriz con los vectores de colocación de los elementos; L vector con las 
longitudes de los elementos; seno, cose, vectores con el valor de las funciones 
seno y coseno de los elementos; nudt es el numero total de elementos; Fm es la 
matriz que contiene el numero de elemento con carga vertical uniforme, el valor de 
la carga, etc.; E es el modulo de elasticidad del material; ELEM matriz con la base 
y la altura de las secciones de todos los elementos; FF matriz con las fuerzas y 
momentos en coordenadas locales de todos los elementos en este orden: fuerza 
axial, corte y momento, del nudo inicial y del nudo final; fact es un factor que ayuda 
a la presentación de los dibujos con los diagramas de momento, corte y carga axial; 
nelem es un vector en el que se indica de que elemento a que elemento se desea 
graficar. 


° dibujo __cortante(X, Y,NI,NJ,q, VC,L,seno,cose,nudt,Fm,E,ELEM,FF,fact,nelem) 


El programa dibujo_cortante encuentra el diagrama de corte; la entrada de 
datos es la misma que la función dibujo_axial. 


e dibujo_momento(X, Y, NI, NJ, q, VC, L, seno, coseno,nudt, Fm,E, ELEM, FF, fact, nelem) 
Con el programa dibujo_momento se halla el diagrama de momentos y los 


datos de entrada son los mismos de la función dibujo_axial. 


e EJEMPLO 2 
Presentar los diagramas de carga axial, corte y momento, utilizando las 


funciones: dibujo_axial; dibujo_cortante; dibujo_momento Estas dos últimas 
funciones grafican los valores de las tablas 2 y 3. 


e SOLUCIÓN 
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Figura 4 a) Diagrama de carga axial, a partir de ordenada de la elastica. 
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Figura 4 b) Diagrama de Momentos; c) Diagrama de corte. 


Desde un punto de vista riguroso, al observar el diagrama de momentos 
mostrado en la figura 4 b, no existe equilibrio de momentos. A manera de ejemplo, 
en el nudo superior izquierdo el momento en la columna es 1.90 Tm; y el momento 
en la viga es 1.83 Tm 


3. DIAGRAMAS A PARTIR DEL EQUILIBRIO 


En el apartado 2 se obtuvo los diagramas de flexión y de corte a partir de 
las ordenadas de la elástica y se vio que desde un punto de vista muy riguroso no 
existe equilibrio en los nudos, falta menos de una décima. Desde el punto de vista 
práctico los resultados obtenidos en el apartado anterior son más que satisfactorios 
más no así desde el punto de vista académico. 


Por este motivo, en el presente apartado se indica como se hallan los 
diagramas, a partir del equilibrio de fuerzas y momentos en los elementos, para ello 
se requiere contar con la solución total de la estructura; en otras palabras, contar 
con las fuerzas y momentos en cada elemento a nivel de los nudos. Para el Ejemplo 
1, que se está analizando en la tabla 4 se indican las cargas resultantes en cada 
elemento. 


Tabla 4 Fuerza axial, de corte y momento en nudo inicial y final de los elementos. 
on |v | uM м мм 

1 

2 

3 
9 


5 





En la tabla 4 se indica la fuerza axial N, el cortante V y el momento M, en el 
nudo inicial y lo propio en el nudo final. Estas cargas, se han colocado en cada uno 
de los elementos, al igual que la carga uniforme distribuida, en la figura 4 d, y se ha 
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procedido a realizar el equilibrio en cada elemento, para encontrar las ecuaciones 
que definen la fuerza axial, fuerza de corte y momento que se muestran en tabla 5. 


2 T/m 2 T/m 


LLTIT ar ia 


0.944 


3.629 


0.944 
= 


0.944 





Figura 5 Fuerzas y Momentos en los nudos de cada elemento y carga uniforme 
distribuida. 


Tabla 5 Ecuaciones de fuerza axial, corte y momento 


M, = 0.989 — 0.944x 


N, = —9.621 
V, = —0.581 


M, = 0.646 — 0.581x 


N, = —4.751 
V, = 1.525 


M, = —1.394 + 1.525х 


М„ = —0.944 
V. = 3.629 — 2x 


2 


2 
М, = —1.843 + 3.629x — > 


М„ = —1.522 


V. = 5.249 — 2х 
2 


2 
М, = —4.427 + 5.249x — D 
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3.1 Nuevas Funciones para CE/NCI-LAB. Segunda Parte 


Para realizar un programa que grafique los diagramas de axial, cortantes y 
momentos para tomando en cuenta el equilibrio de momentos y las fuerzas ya 
calculadas, se considera una ecuacion general, a partir del equilibrio de las cargas 
que actuan en la figura 4 e. 





A 

e | | d y 

ЕАС <<<  _ —_—— = 

мА, MAy 
VI j V' 
(e) | 
Figura 6 Fuerzas y Momentos finales en los nudos y carga uniforme distribuida en 
el elemento. 


Al hacer un cuerpo libre a una distancia x, del nudo inicial y del respectivo 
equilibrio se obtienen las ecuaciones para calcular la fuerza axial N(x), la fuerza de 
corte V (x) y la ecuación de momentos M(x), que se indican a continuación. 

N(x) =—N 
V(x) =V—-q*x 
х? 


M(x) =—-M+V«x—q*—> 


Con estas ecuaciones se realizan los programas para obtener el diagrama 
de axial, cortante y momento, son: 


e dibujo_axial_eq(X, Y,NI,NJ,VC,L,seno,coseno,nudt,Fm,FF,fact,nelem) 
e dibujo_cortante_eq(X, Y,NI,NJ, VC,L,seno,coseno,nudt,Fm,FF,fact,nelem) 
e dibujo_momento_eq(X, Y,NI,NJ,VC,L,seno,coseno,nudt,Fm,FF,fact,nelem) 


Los datos para estos programas son los mismos que los indicados en el 
apartado 2.4, por lo que se omite esa informacion. 


En las figuras 4 f, g, h, se indican los diagramas de fuerza axial, de corte y 
de momentos, que se hallan a partir del equilibrio de los elementos. 
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Diagrama de Cargas Axiales 





(a) 


Diagrama de Corte 
5.25 





(b) 


Diagrama de Momento 





(c) 
Figura 7 Diagramas obtenidos a partir del equilibrio de elementos; a) Axial; b) Corte; 
c) Momentos. 
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4. FUNCIONES DE FORMA PARA CARTELA A LA DERECHA 


En Aguiar y Bonilla (2000, 3) se presento el calculo de funciones de forma 
que conducen al calculo de las acciones de empotramiento perfecto, para el 
elemento cartela a la izquierda, no se lo realizó para el elemento cartela a la derecha 
debido a que se está trabajando el caso de vigas con dos cartelas simétricas. Ahora 
que se va a dibujar el diagrama de momentos y corte, es necesario tener definidas 
las funciones de forma para el elemento indicado en la figura 5, que corresponde a 
elemento cartela a la derecha. Nótese que se está trabaja con Modelo 3, en donde 
solo se requiere la variable a. 





X 
Figura 8 Elemento Cartela a la derecha. 


Las ecuaciones para este caso de cartela se realizan de la misma manera 
que en caso de cartela izquierda, comenzando con una carga Po=0, la cual se 
utilizará para encontrar la ecuación de momentos de la que parte la función de 
forma. 


OMS р E 
d? x = 
ам 
ле (16) 
ах н 
М(х) =А*х+В (17) 


Para encontrar la ecuación de la elástica v(x) se parte de la siguiente 
ecuación diferencial, teniendo en cuenta que la variación del momento de inercia es 
variable. 


div М(х) 
dx?" EI(x) 





La ecuación de la variación de la altura de la cartela (ver figura 5) viene 
dada, por: 


H-h 


x=0;y=h == 
x=Ly=H у= ћ = m(x — 0) 
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H—h 





y=h+ 
De lo que la inercia de la viga es: 


¡ee — | 
190 p 


X 








— 
Хх 


р 
I) => (aH + Т 


З З 


| _b*H ( Le 
о" a+t——*x 








L 
РН? 
mia 
1-а 1 
I(x) = Í. (a+ L ex) (18) 


Reemplazando en la ecuación de resistencia de materiales se tiene: 


d?v А*х+В 


ЕУ; = 3 


E * la * (a += х) 





Resolviendo la ecuación diferencial utilizando el método de integrales 
sucesivas se tiene: 


dv L?((2Ax — AL + В)а – 2Ах – В) 


ах 2E l, (æ — 1) (ха – Іа — х)? 9) 


2AIn(x(a —1)—La) (AL + В)а – В ) 
а – 1 x(a? — 24 + 1) - La? + La 
2Е1, (а – 1)? 


1? ( 
v(x) = + Cx + D 
(20) 
4.1 Función de forma %0;(x) 


Solo existe v, en el nudo inicial. Con esto las condiciones de borde son: 


I) x=0 0=0 
li) x=0  v(x)=v1 
lil) х= 0=0 
iv) x=L v=0 


Con estas condiciones de borde, las constantes de integración son: 
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El, v,a*—2E1, v,a7+2E1,v,a—-El,v, 


ar 13 a@(2 — In(@)) + L3(—In(a@) — 2) 


В El, v, (a— Da 
~ L2((a + 1) In(a@) — 2a + 2) 
сс v, (a— 1) 
~ L(aln(a) + In(a) — 2a + 2) 
— ,(In(—L)a + 2a + In(-L) — 1) 
аа) + (а) — 2a + 2 


Reemplazando estas constantes en le ecuación 20, se obtiene: 


0,(х) 
_ (Pa? + La + Lx) In(- a =) + Lxa? eee 


=D> ) + (xx? — 3Lx + 2L2)a? + (-2x? + 4Lx — 2L5a + x?— Lx 
L((x —L)a—x)((a+1)In(@) — 2a + 2) 


(21) 


42 Función de forma 0x(x) 


Solo existe 6, en el nudo inicial. Las condiciones de borde son: 


I) x=0 Ө = 01 
li) x=0 v=0 
111) v= 0=0 
iv) x=L v=0 


Resolviendo el sistema de ecuaciones se obtienen las constantes de 
integracion: 
vee (El, at-3El,04+3E1,0?—E1,0 0, 
o L2(In(a) + 2) + L2a(In(a) — 2) 
ЕТ. а? (210(а) + а? – 40 + 3) 0, 
o L(aln(a) + In(a) — 2a +2) 
E (a7(In(a) — 1) +a) 0, 

 —In(oa) + a2(In(a) — 2) + 4а – 2 


L a (a In(L7a@7) — a7 In(—La) — a — In(—L) + 1)0, 


= (а — 1)2((« + 1)ш(о) — 2« +2) 


Al sustituir las constantes de integracion en la ecuacion 20, se tiene: 


Ø (x) 
a aS + a 9%) 4 py DO — y2@ n(a) + (x? — Lx) — x? + Lx) 
((x — L)a — x)((a + 1) In(@) — 2a + 2) 


(22) 
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4.3 Calculo de @-(x) 


Solo existe v,; las condiciones de borde son: 


I) x=0 0=0 
11) x=0 v=0 
111) x= L 0=0 
iv) x=L v=v2 


Al resolver el sistema de ecuaciones, que se forman con las condiciones 
de borde, se encuentran las constantes de integracion. 


E l, v (a — 1)?’ (a + 1) 


AS e aa 


ee E l, v, (a — 1)°a 
~ 12 (a In(a) + In(a) — 2a + 2) 


v (æ — 1) 


nye na) 0) 


_ v2 (aln(—La) + In(—La) + 1) 
~~ aln(a) + In(a) — 2a +2 


Al reemplazar las constantes de integración en la ecuación 20, se halla. 


= De — (x? — Lx)æ? + 2x?’ a — x? — Lx 


(La + La + Lo n= Ga ~ txa? n@— =e 
Ds(x) = L((x — L)a — x)((a + 1) In(a@) — 2a +2) 
(23) 


4.4  Calculo de @,(x) 


Solo existe el giro 8, en el nudo final. Con esto las condiciones de borde 
son las siguientes. 


I) x=0 0=0 
11) x=0 v=0 
111) Soh 0 = 02 
iv) x=L v=0 


Al resolver el sistema de ecuaciones, que se forman con las condiciones 
de borde se obtiene: 


Ë To (z — 1)? 0, 


АС Соо) 


_ E L, a(2a In(a) — a* +1) 0, 
~ L(aln(a) + In(a) — 2a + 2) 
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_ (—In(a) + a — 1) 6, 
(а) + a?(In(a) — 2) + 4a — 2 


L(a In(LZa) — æ? In(—La) — In(—La) + a — 1) 0, 
o (a — 1)2(aIn(a) + In(a) — 2a + 2) 
Al sustituir las constantes de integración en la ecuacion 20 se halla 0, (x) 


Q G) 


(La + Lx) Int G=Da=% 


((х = Ша – х)((а + то — 2a + 2) 


а L)a 


) + Lxa nt ) — x? In(a) + (x? — Lx)a — x? + Lx 


(24) 


5. FUNCIONES PARA DIBUJAR DIAGRAMAS 


Para encontrar los diagramas de momento, corte y carga axial, se debe 
determinar primero las ordenadas de la elástica u(x),v(x) teniendo presente que 
estas tienen dos contribuciones a saber: la primera del problema complementario 
(funciones de forma) y la segunda del problema primario que consiste en resolver 
la viga empotrada-empotrada ante cualquier tipo de carga; en este artículo se 
desarrolla el caso de carga uniforme distribuida para los dos casos: cartela a la 
izquierda y cartela a la derecha; el desarrollo numérico se presenta en el Anexo 1. 


Para la deducción se parte de la ecuación 10, la misma que se repite a 
continuación. 


V(X) = 0,P2(x) + 0, Pz (x) + „ф»(х) + 0,ó, (x) + Vy 


Para vigas de sección constante se tienen las funciones de forma indicadas 
en el apartado 2 y para vigas con cartelas (para el modelo 3) las funciones de forma 
se encuentran en el apartado 4, para el caso de cartela a la derecha ya que para 
cartela a la izquierda se dedujeron en Aguiar y Bonilla (2020, 3). Es importante 
destacar que esto representa el Problema Complementario. Para el Problema 
Primario se va a indicar como se lo hace mas adelante y en el Anexo 1. 


5.1 Diagrama de Corte 





dx? dx? dx? dx? dx? dx 


аЗъ(х) = d° (ui 2 (x)) + а2(0,фз(х)) + d(v,ps(0)) + 43(Ө›фь (х)) _ dš (w) р) 
Se ha colocado el signo negativo en el ultimo termino, por la convencion de 
signos, que ya se indicó anteriormente. 


Al tener en cuenta que los desplazamientos v,,v, y los giros 6,,0., son 
constantes, se tiene. 
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divo аф), а(@3(0)) P), n E) dlo) 
dx? ~ "4 dx? T dx? e dx? a ахз ахз 





Entonces se debe resolver la viga empotrada-empotrada соп carga 
uniforme distribuida y hallar v, para posteriormente derivar tres veces y encontrar el 
cortante del Problema Primario. De esta forma no se lo hace porque es demasiado 
largo se llega al cortante recordando que la derivada del corte V es igual a la carga 
vertical P, 


dv 
dx “2 
Al integrar se tiene: 
V= P, *x+A (25) 


Al integrar la ecuación 25, se halla el momento (no se debe perder de vista 
que se está resolviendo el Problema Primario para el caso de viga empotrada- 
empotrada con carga uniforme P,) 


2 
Po x 





M = +Ax+B (26) 


El cálculo de las constantes de integración A, B, se indica en el Anexo 1. Se 
destaca una vez mas que la ecuacion (25) reporta el cortante del Problema Primario 
y la ecuacion (26) el momento del Problema Primario. A continuacion, se presentan 
estas constantes en forma ejecutiva, para el Modelo 3. 


e Para cartelas izquierdas: 


_ L Po (4a In(@) + 2 In(@) – а? – 4а + 5) 


(27) 
2(a —1)(a In(a) + In(a) — 2a + 2) 
e Para cartelas derechas: 
L P, a(—In(a*) — 4) + L P, œ? (5 — In(a?))— LP 
is ee 0 ( ) 0 (28) 


In(a*) + a (4 — log(a?)) — 8a +4 
Para el caso de la sección constante de la viga es mejor derivar tres veces 
la ecuación 12, que se repite a continuación (Ver apartado 2.3) 


Po 
24El 





V(x) = 0, P2(x) + 6,3 (x) + >›ф»5(х) + 0;ф„(х) — (X* + Х^1° — 2ХЗГ) 
5.2 Diagrama de Momentos 


Para la construcción de diagramas de momento es necesaria aplicar la 
misma lógica anterior solo que para la segunda derivada. 
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v(x) = vıp: (x) + 0,3 (x) + Vzbs (x) + O26 (x) — Vy 


Pre) _ Ф (фо), 4 (0,фа(х)), а (ооф), @(0әфьб)) @(%) 


ах? ах? ах? ах? ах? ах? 


Al tener en cuenta que los desplazamientos v,,v, y los giros 0,,0,, son 
constantes, se tiene. 


ano 2. Plat) PO) LP) Po) 
dx? TA dx? т dx? e dx? +8, dx? ах? 





2 2 


Vy x 
dx? *EI =M = Pox; +A*x+B 





Para tener completo el formulario se vuelve a copiar la constante de 
integración A, para los dos casos. 


e Para cartelas izquierdas: 


_ LP, (4a In(a) +2 In(a) – а? – 4а + 5) 
-= 2(a —1)(a In(a) + In(a) — 2a + 2) 


B 
_ 2 Po (In (L) о? In (1? <“) +a? In(L? cŠ) — 2 x? In(L <)2 + (—2 « – 1) т(о) +07 + (—2 In(L)* — 3 In(L) — 5) «?+ 7 < —3) 
р 2(x< —1)3((« +1) In(x) — 2 « +2) 


(29) 


e Para cartelas derechas: 


_LP, a(—In(a*) — 4) + LP, a?(5 — In(a?)) — LP, 


A= 
Inla?) + a? (4 — log(a?)) —8a+4 


B 
_ L?Py & (In(—L) « In(L* &*) —0 In(—L? š) — 2 < In(—L x)? + (x?+ 2 x?) In(x) — 3 «P+ 7 «+ (—2 In(—L)* + 3 In(-L) — 5) « +1) 
o 2(« —1)3((« +1) In(«) — 2 « +2) 


(30) 


6. NUEVAS FUNCIONES DE CEINCI-LAB. PARTE 3 


La entrada de datos es muy similar a la indicada en el apartado 2.4, por lo 
que a continuación se presentan los programas que obtienen los diagramas de 
carga axial, corte y momento; y a continuación se describe brevemente como lo 
hace. 
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Pero hay una parte comun para los tres diagramas, cual es encontrar para 
cada elemento los valores de las coordenadas locales u,, v,,6,,u,, v,,0,. Esto se lo 
realiza de la siguiente manera: 


1) Con el vector q y con la matriz que contiene a los vectores de colocación 
VC se determina el vector de deformaciones en coordenadas globales p, 


li) Se determina la matriz de rotación T», y se encuentra el vector de 
deformación en coordenadas locales p, 


Pi = 123 * Dg (31) 


lil) Los elementos del vector p, son los valores de uj, v4, 04, Us, V2, 0, 


e dibujo_axial_cartelam3(X, Y,NI,NJ,q, VC,L,seno,coseno,nudt,Fm,E,ELEM, 
FF, fact,nelem) 


Una vez que se tienen los valores de u, v4, 01, u2, V2, 02, en cada elemento, 
se encuentra u(x) trabajando con la ecuación 9, previamente se define que x es 
una variable simbólica (syms x). 


Se deriva u(x) con la función de Matlab diff. Luego se obtiene el valor de la 
fuerza axial Nen 10 puntos del elemento. El programa determina el área transversal 
de la sección a partir de los dos primeros valores del arreglo ELEM; cuando se trata 
de un elemento circular el área es el primer valor de ELEM 


e  dibujo_momento_cartelam3(X, Y, NI, NJ,q,VC, L, seno,coseno,nudt, Fm, E, EL 
EM, FF, fact, nelem) 


El archivo ELEM contiene b, h, H, tipo; que son el ancho de la viga, la altura 
al inicio de la cartela, la altura al fin de la cartela y la variable identificadora que vale 
0, si cartela se encuentra a la izquierda, 1, si esta a la derecha y 2 si es de sección 
constante. 


En base a esta variable identificadora se encuentran las funciones de forma, 
para cartela a la izquierda las funciones de forma se encuentran en Aguiar y Bonilla 
(2020, 3); para cartela a la derecha las funciones de forma son las deducidas en 
este artículo y para sección constante las indicadas en apartado (2.1). 


Para el caso de cartela se tiene: 





v(x) = 0,P2(x) + 01фз(х) + у,ф5(х) + 0,ф,(х) – = + Ах + в) 


Si la cartela se halla a la izquierda se trabaja con las ecuaciones 27 у 29; 
para el caso de cartela a la derecha con ecuaciones 28 y 30. 
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e dibujo_cortante_cartelam3(xX, Y,NI,NJ,q, VC,L,seno,coseno,nudt,Fm,E,EL 
EM, FF, fact, nelem) 


e EJEMPLO 3 


La estructura que se analiza tiene columnas cuadradas de 1.20 m, de lado 
y vigas con cartelas lineales simétricas de 0.5/0.7 m, en el tramo de sección 
constante y 0.5/1.0 m, en el extremo de la cartela, ver figura 5 a; la altura del primer 
piso es 4.5 m, y de los pisos superiores es 4.0 m. En el primer y tercer vano de 8.0 
m, la carga uniforme es de 2.80 T/m, igual en todos los pisos y en el vano central 
de 10.0 m, la carga es de 3.45 T/m; ver figura 5 b. El Módulo de elasticidad del 
hormigón es E,, = 2400000 T /m?. Se pide realizar el análisis estático con el modelo 
3, para lo cual en la figura 5 c, se indica la numeración de los nudos y elementos. 
Se pide presentar los diagramas de momento, corte y carga axial, empleando las 
funciones de CEINCI-LAB y para el primer piso indicar las fuerzas y momentos en 
los elementos. 


W4= 2.80 T/m W4= 3.45 T/m W4= 2.80 T/m 
PT LILAC eee Ee LCL. Lo 








1.20/1.20 1.20/1.20 


к W3= 280 T/m| W3= 3.45 T/m W3= 2.80 T/m 


W2= 2.80 T/m W2= 3.45 T/m W2= 280 T/m 


1.20/1 20 1.20/1.20 
0.70 





4.00m .4.00m . 4.00m 


Wi= 280 T/m Wi= 3.45 T/m Wi= 2.80 T/m 


1.20/1.20 1.20/1.20 
0.70 


1.20/1.20 1.20/1.20 


4.50m 
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° 1 °2 °з @4 
(с) 
Figura 9 Estructura de 4 pisos con columnas cuadradas y vigas con cartelas de 
hormigón armado; a) Geometría del Pórtico; b) Cargas verticales; c) 
Numeración de nudos y elementos. 
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e SOLUCIÓN 


El programa desarrollado se denomina: Ej 77 Estructura 4 pisos se 
encuentra en la carpeta cartelas en la siguiente dirección electrónica: 
https://www.researchgate.net/profile/Roberto Aguiar2/research 


o) 
1509 


16485 





(b) 


Figura 10 a) Diagrama de Momentos; b) Diagrama de Corte 
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Cargas Axiales 





) 
Figura 10 c) Diagrama de Carga Axial 


En el calculo se obtiene en primer lugar la contribucion a la matriz de rigidez 
de las columnas cuadradas, empleando para el efecto el programa krigidez, que a 
su vez llama al programa kmiembro para encontrar la matriz de rigidez de los 
elementos en coordenadas globales. 


Después para las vigas se obtiene la contribución a la matriz de rigidez con 
el programa krigíidez_cartela_dosnudosintermedios por ensamblaje directo; este 
programa llama a: kmiembro_cartela_dosnudos intermedios para hallar la matriz 
de rigidez en coordenadas locales. Este último programa mediante la variable 
indicadora tipo obtiene la matriz de rigidez de la cartela izquierda si esta variable 
indicadora vale 0; de la cartela derecha si vale 1; de la parte intermedia de sección 
constante si vale 2. Estos datos al igual que la geometría de la cartela se dan en el 
arreglo ELEM. 


Una vez que se obtienen los desplazamientos y giros de la estructura, se 
halla las fuerzas en las columnas con el programa fuerzas y en las cartelas con el 
programa fuerzas_cartela_modelo_3. Con estos programas se determina la fuerza 
axial, de corte y momentos en el nudo inicial y final de los elementos en coordenadas 
locales; los mismos que se agrupan en una sola matriz para dibujar los diagramas 
de momento (figura 6 a) de corte (figura 6 b) y de carga axial (figura 6 c). En Aguiar 
y Bonilla (2020, 1, 2, 3) se ha indicado a mas de los modelos de analisis de vigas 
con cartelas, se han indicado una gran cantidad de funciones de CEINCI-LAB. De 
igual forma hay información en Aguiar (2020) y en el blog de CEINCI-LAB. 


https://ceincilab.wordpress.com 


En la figura 7 a, se presentan las fuerzas y momentos en cada uno de los 
elementos del vano izquierdo de la planta baja; en la figura 7 b, se tiene del vano 
central y en la figura 7 c, del vano derecho. Al ser la estructura y la carga simétrica 
los resultados también deben guardar esta simetría. Estos resultados reportan los 
programas fuerzas y fuerzas_cartela_modelo_3, en una matriz; el usuario debe 
saber interpretar y obtener los resultados en cualquier elemento. 
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(c) 
Figura 11 Fuerzas y momentos en los elementos de la planta baja; a) Vano 
izquierdo; b) Vano central; c) Vano derecho. 
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Con un poco de detenimiento se pueden ver estos resultados en los 
diagramas de momento, corte, indicados en la figura 6. 


e VERIFICACIÓN CON SAP2000 


Los resultados encontrados con el sistema de computación CEINCI-LAB 
fueron comparados con los que reporta el programa SAP 2000. 





(a) 





| | | 
I h зт х ap akh 


(b) 
Figura 12 Comparación de resultados con SAP 2000; a) Geometría que reporta el 
programa; b) Diagrama de Momentos. 
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(d) 
Figura 12 Diagramas obtenidos con SAP 2000; c) Diagrama de Corte; d) Diagrama 
de Carga Axial. 


Existe muy buena relación entre los diagramas de flexión, corte y carga axial 
que reporta CEINCI-LAB con los diagramas encontrados con SAP 2000; esto valida 
los modelos que se han desarrollado sobre cartelas de hormigón armado. 
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e EJEMPLO 4 


La estructura de la figura 9 a, tiene columnas circulares de 1.20 m, de 
diámetro y las vigas rectangulares son de 0.5/0.7 m; en el vano central que es de 
11 m, se tiene una cartela de 2.0 m, en los extremos y la altura en el extremo es de 
1.10 m. En la figura 9 b, se tienen las cargas verticales que gravitan en la estructura 
en los vanos de 5.5 m, la carga vale 3.89 T/m, del piso1 al 4, en el piso 5 la carga 
es de 3.85 T/m; en el vano central de 11.0 m, la carga vale 6.53 T/m, del piso 1 al 4 
y en el piso 5 es de 6.45 T/m; el módulo de elasticidad E, = 2400000 T /m?. Se pide 
presentar los diagramas de momento, corte y carga axial con las funciones de 
CEINCI-LAB y comparar los resultados con el programa ETABS 
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Figura 13 Estructura con vigas con cartelas de hormigon armado; a) Dimensiones 
de las vigas y columnas; b) Carga vertical actuante; c) Numeracion de 
nudos y elementos. 
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e SOLUCION 


En la figura 9 c, se indica la numeracion de nudos y elementos con que se 
resuelve el Ejemplo. Notese que en los vanos exteriores que no se tienen cartelas 
se ha considerado un solo nudo en la mitad del elemento y en los vanos con cartelas 
dos nudos para trabajar con el modelo de cartelas número 3. 

El programa que resuelve el Ejemplo es: Ej_10_Estructura_5_ pisos 


Para el caso de columnas circulares la contribución a la matriz de rigidez de 
la estructura se realiza con el programa krigídez_acero ya que se dan como datos 
el área y la inercia de la sección transversal. Por otra parte, para encontrar las 
fuerzas y momentos en los elementos se debe utilizar el programa fuerzas_acero 
por el mismo motivo. 


Para obtener los gráficos de los diagramas de momentos, corte y carga axial 
se debe trabajar con una matriz ELEM diferente a la que se consideró la contribución 
a la matriz de rigidez. Ahora el contiene el radio de la columna circular, se vuelve a 
repetir está cantidad, después el código para columnas circulares es 3, por lo que 
se da dos veces ese dato; en total se tiene 4 datos por cada elemento, que es igual 
al número de datos de las vigas con cartelas. En resumen, para hallar los diagramas 
de momento, corte y carga axial, el código es 2 cuando las columnas son 
rectangulares o cuadradas y 3 cuando son circulares. 


Momentos 





0 


10 19 20 


(a) 
Figura 14 a) Diagrama de Momentos con funciones de CEINCI-LAB 


632 Roberto Aguiar, Alejandro Bonilla y Jorge Romero 





et | a i — [лын шы F п ad 
=т=т  — z UCE y = i 


(b) 
Cargas Axiales 
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(c) 
Figura 14 b) Diagrama de corte; c) Diagrama de carga axial, con funciones de 
CEINCI-LAB 
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Figura 15 a) Diagrama de Momentos 
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(b) 





(с) 
Figura 15 b) Diagrama de Corte. Encontrados con ETABS. c) Diagrama de carga 
axial obtenido con programa ETABS. 
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ANEXO 1 


PROBLEMA PRIMARIO EN VIGAS CON CARTELAS 


Para el problema primario en cartelas se debe considerar la ubicacion de la 
cartela (si esta a la izquierda o si esta a la derecha) y colocarla como biempotrada 
con una carga uniformemente distribuida P}. 


Para llegar a la ecuacion de deflexiones, es necesario llegar a la ecuación 
de momentos: 
d? M(x) 
ах? ° 
Por lo que se debe recurrir al método de integrales sucesivas para resolver 
esta ecuación: 
dM (x) 
dx ` 





Po *x+A (A.1) 


2 
x 
М(х) = + +А*х+В (А.2) 
Se pide al lector que imaginariamente piense en una viga sometida a una 
carga uniforme P, y en el nudo inicial se tiene un cortante de magnitud A, y un 
momento de magnitud B. La ecuación A.1 se halla por equilibrio de fuerzas 
verticales y la ecuación A.2 por equilibrio de momentos, a una distancia x. 


A partir de este punto, se incluye el efecto del material y la sección en el 
cálculo, para esto, se recurre a la ecuación de resistencia de materiales. 





div М(х) 
dx2  El(x) (A.3) 
e Para cartela izquierda 
х? 
2 Po x= +A +B 
a u = s (A.4) 
d = 
di Exlax(1- ax) 


Recurriendo nuevamente al método de integrales sucesivas, se tiene: 


8(x)=(L43*((2*Po*x*2*a42+(4*L*Po*x-4*Po*x2)*q+2*Po*xh2- 
4*L*Po*x+2*L42*Po)*In(x*a-x+L)+(-4*A*x-2*B)*a%2+((4*L*Po+8*A)*x- 
2*A*L+4*B)*a+(-4*L*Po-4*A)*x+3*L*2*P0+2*A*L-2*B))/(4*E*la*(a-1)*3*(x*a- 
X+L)42)+C 


(A.5) 
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v(x)=(((2*LA3*Po*xA2-4*A*LA3*x)*qA2+(-4*LA3*Po*xA2+(8*LA4*Po+8*A*LA3)*x- 
4*A*LA4)*o+2*LA3*Po*xA2+(-8*LA4*Po-4*A*LA3)*x+6*LA5*Po+4*A*LA4)*In(x*o- 
x+L)+4*C*E*|a*xA2*oA5+(4*C*E*|a*L*x-20*C*E*la*xA2)*oA4+(40*C*E*la*xA2- 
16*C*E*la*L*x)*oA3+((-2*LA3*Po- 
40*C*E*|a)*xA2+24*C*E*|a*L*x+2*B*LA3)*oA2+((4*LA3*Po+20*C*E*la)*xA2+(- 
2*LA4*Po-16*C*E*la*L)*x-2*A*LA4-4*B*LA3)*o+(-2*LA3*Po- 
4*C*E*la)*xA2+(2*LA4*Po+4*C*E*la*L)*x+LA5*Po+2*A*LA4+2*B*LA3Y/(4*E*la*(o- 
1)^4*(x*a-x+L))+D 

(A.6) 


Se tienen 4 constantes de integración: A, B, C, D, y 4 condiciones de 
frontera: 


i) x=0 0=0 
11) x=0 v=0 
111) x=L @=0 
iv) x=L v=0 


Aplicando estas condiciones, se resuelve las constantes de integracion teniendo que 
estas son: 


A=(Po*(L*(-2*In(L*a)-5)+2*L*In(L))+Po*a*(L*(4- 
4*In(L*a))+4*L*In(L))+L*Po*a"2)/(2*In(L*a)j+ar2*(-2*In(L*a)+2*In(L)+4)-8%*a- 


2*In(L)+4) (А.7) 


В=-(Ро*а^2*(1^2*(-2*|п(1*в)^2+3*1п(1*ов)-5)+^2*|п(1)*(4*1п(1*о)-3)- 
2*LA2*|n(L)42)+Po0*(L*2*(-In(L*a)-3)+L42*In(L))+Po* a*(L42*(7- 
2*In(L*a))+2*L12*In(L))+L2*Po*a13)/(a73*(4*In(L*a)-4*In(L)- 
16)+2*In(L*a)+a14*(-2*In(L*a)+2*In(L)+4)+a*(-4*In(L*a)+4*In(L)-16)+24*a12- 


2*In(L)+4) [д] 


C=-(Po*oA2*(LA3*(-In(L*a)A2+2*|n(L*o)- 
1)+LA3*|n(L)*In(L*oa))+Po*(LA3*|n(L)*(In(L*a)+2)-LA3*|n(L)A2-LA3)+Po*o*(LA3*(2- 
2*In(L*a))-2*L3*In(L)))/(a*(E*la*(6*In(L*a)+20)- 
6*E*la*In(L))+a14*(E*la*(6*In(L*a)-20)-6*E*la*In(L))+E*la*(-2*In(L*a)- 
4)+a%5*(E*la*(4-2*In(L*a))+2*E*la*In(L))+a%3*(E*la*(40- 
4*|n(L*a))+4*E*la*ln(L))+a%2*(E*la*(-4*In(L*a)-40)+4*E*la*ln(L))+2*E*la*In(L)) 


D=(Po*a%2*(L44*(-In(L*a)*%2-2)+L%4*In(L)*(In(L*a)- 
3))+Po*o*(LA4*|n(L)*(2*In(L*a)+3)-2*LA4*|n(L)A2+LA4)+Po*(LA4*|n(L)*(-In(L*o)- 
1)+104*In(1)42)+(144*In(L)+L14)*Po*a73)/(a12*(E*la*(10*In(L*a)+60)- 
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10*E*la*In(L))+ar5*(E*la*(8*In(L*a)-24)- 
8*E*la*In(L))+E*la*(2*In(L*a)+4)+a"6*(E*la*(4- 
2*In(L*a))+2*E*la*In(L))+a*(E*la*(-8*In(L*a)-24)+8*E*la*In(L))+an4*(E*la*(60- 
10*In(L*a))+10*E*la*In(L))-80*E*la*a”3-2*E*la*In(L)) 

Se ha numerado solo las ecuaciones de las constantes A, B porque esas 
son las que interesan para encontrar el cortante y momento con las ecuaciones 
A.1 y A.2 


° Para cartela derecha 





2 

421 Pr TA*x+P 
dx2 f= 9 
E + la» (a+ + x) 


Al integrar se encuentra el giro y el desplazamiento vertical 


8(x)=((x*(4*A*L43*042+(4*L44*P, *a-8*A*L43)*a- 

4*L*4* Py *at4*A*L43)+2*B*L43*0%2+(-2*A*L*4*a-4*B*L43)*a- 

З^“| ^Б*Р„*а^2+2*Д*| ^4*а+2*В* ^3З)/(х^2*(4*а^5-20*а^4+40*а^3-40*а^2+20*а- 
4)+х*(-8* *а*0^4+32*| *а*0^3-48*| *а*0^2+32*| *а*а-8**а)+4*|^2*а^2*а^3- 
12*LA2* ar2*ah24+12*LA2% ah2*a-4*L*2*a%2)-(L*3* P, “In(x*(a-1)-L*a))/(2*a43- 
6*a*2+6*a-2))/(E*la)+C 


v(x)=(-(2*B*L^3*a^2+(2*A*L^4*a-4*B*L^3)*a+L^5*Po*a^2- 
2*A*L^4*a+2*B*L^3)/(x*(4*a^5-20*a^4+40*a^3-40*a^2+20*a-4)- 
4*L*a*a^4+16*L*a*a^3-24*L*a*a^2+16*L*a*a-4*L*a)+(In(x*(a-1)- 
L*a)*(A*L^3*a+L^4*P,*a-A*L^3))/(a^4-4*a^3+6*a^2-4*a+1)-(L^3*P,*(-x*(0- 
1)+In(x*(a-1)-L*a)*(x*(a-1)-L*a)+L*a))/((a-1)*(2*a^3-6*a^2+6*a-2)))/(E*la)+C*x+D 


Las condiciones de borde son: 


I) x=0 0=0 
11) x=0 v= 
lil) x= L Ө = 
iv) x=L v=0 


Aplicando estas condiciones, se resuelve las constantes de integración teniendo 
que estas son: 
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A=-(L*Po*a42*(-2*|n(L*(a-1)-L*a)-5)+L*Po*a*(4-4*|n(L*(a-1)- 
L*a))+(2*L*Po*a42+4*L*Po*a)*In(-L*a)+L*Po)/(2*In(L*(a-1)-L* a)+a%2* (- 


2*|n(L*(a-1)-L*a)-4)+(2*a42-2)*In(-L*a)+8*a-4) (A.9) 


B=-(L42*Po*a*2*(2*In(L*(a-1)-L*a)*2-3*In(L*(a-1)- 

L* a)+5)+L42*Po*a*3*(2*In(L*(a-1)-L* a)-7)+L*2*Po*a*4*(In(L*(a-1)-L* a)+3)+In(- 
L* a)*(LA2*Po*ah2*(3-4*In(L*(a-1)-L* a))-L42*Po*a%4- 
2*LA2*Po*ah3)+2*LA2*Po*ah2*In(-L* a)42-L42*Po* a)/(a*3*(4*In(L*(a-1)- 
L*a)+16)+2*In(L*(a-1)-L*a)+04*(-2*In(L*(a-1)-L%a)-4)+a*(16-4*In(L*(a-1)- 


L*aq))+(2*a14-4*013+4*a-2)*In(-L*a)-24*a1 2-4) 
(A.10) 


C=(L43*Po*(In(L*(a-1)-L*a)2-2*In(L*(a-1)-L*a)+1)+In(-L*a)*(-L43*Po*In(L*(a-1)- 
L*a)+113*Po*a"2*(-In(L*(a-1)-L%a)-2)+2*L43*Po*a)+L3*Po*a*(2*In(L*(a-1)- 
L%a)-2)+L3*Po*a2*In(-L*a)2+L43*Po*a12)/(E*la*ar4*(6*In(L*(a-1)- 
L*a)+20)+E*la*a*(6*In(L*(a-1)-L*a)-20)+E*la*(4-2*In(L*(a-1)-L*a))+E*latars5*(- 
2*In(L*(a-1)-L*a)-4)+E*la*ar2*(40-4*In(L*(a-1)-L*a))+E*la*ar3*(-4*In(L*(a-1)- 
L*a)-40)+(2*E*la*ar5-6*E*la*tan4+4*E*la*an3+4*E*la*a”2- 
6*E*la*a+2*E*la)*In(-L*a)) 


D=-(L44* Po* a*(-In(L*(a-1)-L* a)*2-In(L*(a-1)-L*a))+In(- 

L* a)*(LA4* Po*ah2*(2*In(L*(a-1)-L* a)+3)+L44*Po*a*(In(L*(a-1)-L*a)- 
2)+L44*Po*a3*(-In(L*(a-1)-L*a)-2)+L*4*Po)+L*4* Po* a3 *(In(L*(a-1)- 

L* q)+2)+(L44* Po*a%3-2*L*4* Po*a2)*In(-L*a)42- 
3*LA4*Po*aQh24+L44*Po)/(E*la*a*2*(10*In(L*(a-1)-L*a)- 
60)+E*la*ar5*(8*In(L*(a-1)-L*a)+24)+E*la*(2*In(L*(a-1)-L*a)-4)+E*la*ar6*(- 
2*In(L*(a-1)-L*a)-4)+E*la*a*(24-8*In(L*(a-1)-L*a))+E*la*ar4*(-10*In(L*(a-1)- 
L*a)-60)+(2*E*la*ar6-8*E*la*ar5+10*E*la*ar4-10*E*la*an2+8*E*la*a- 
2*E*la)*In(-L*a)+80*E*la*a”3) 


Las constantes arriba presentadas se encuentran en el apartado 5, de una 
forma simplificada para facilidad de lectura. 


